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Abstract 

This work is concerned with existence of weak solutions to discon¬ 
tinuous stochastic differential equations driven by multiplicative Gaus¬ 
sian noise and sliding mode control dynamics generated by stochastic 
differential equations with variable structure, that is with jump nonlin¬ 
earity. The treatment covers the finite dimensional stochastic systems 
and the stochastic diffusion equation with multiplicative noise. 


1 Introduction 

We consider here stochastic differential equations of the form 

dX + AXdt + f{X)dt = B{X)dW, te(0,T) 

V(0)=a:, 

where A: D{A) G H ^ H is self-adjoint, positive definite such that 

is of trace class for some 5 G (0,1), lU is a cylindrical Wiener process of the 

form 

OO 

W{t) = J2nmei- ( 1 - 2 ) 

J=1 

Here {ej} is an orthonormal basis in H, Acj = Xj ej and is a mutually 

independent system of Brownian motions in a probability space {12, P} 
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with filtration The operator fH H is Borel measurable and 

locally bounded while B G L{H, C‘^{H)) where C?{H) is the space of Hilbert- 
Schmidt operators on H. 

It should be said that under these general conditions equation (1.1) is 
not well posed except the case of additive noise {B{X) = I) where (1.1) has 
a unique weak (martingale) solution, see [5]. Equation (1.1) has however a 
unique strong solution if / is Lipschitz or accretive and continuous, see [4], 
or more generally if f is a maximal monotone graph in R x R with domain 
D{f) = R (see [2]). 

Equations of the form (1.1) with discontinuous / describe systems with 
variable structure and, in particular, closed-loop control systems with “slid¬ 
ing” mode behaviour. Here we shall study from this perspective two special 
cases. 

The first one is the finite dimensional system 

dX + f{X) dt = a{X) dW 

X{0) = x ^ 


where W is a. n-dimensional Wiener process and / G Lg,(E",E"), 6 

Lip{W, /.(R"-, R"-)). 

The second one is the stochastic partial differential equation 

dX - AX dt + f{X) dt = b{X) dW, in (0, T) x O 
X = 0, on{0,T)xdO (1.4) 

W(0,O = a;(O, eeO 

in a bounded and open domain O C R'^, d > 1 with smooth boundary dO. 
Here IE is a cylindrical Wiener process of the form (1.2) in H = L^{0), where 
{Pj}’^i is a system of independent Brownian motions in a hltered probability 
space (H, X, Xt, P), fij G R, j = 1, 2,... and {cj} is an orthonormal basis in 
L‘^{0) to be made precise later on. Here / G and b G Lip,„^(R). 

Like in deterministic case, for existence in equation (1.4) one must extend 
it to a multivalued stochastic equation of the form 

dX -AXdt + F{X) dt 3 b{X) dW, in (0, T) x O 
X = 0, on{0,T)xdO (1.5) 

X(0,O = a;(O, eeO 


where F: R —)■ 2^ is the Filippov map associated with /, that is (see [6], [7]) 


F{r) = [m{fr),M{fr)], Vr G R 
m{fr) = lim ess inf f{u) 

(5—>-0 4iG[r—(5,r+(5] 

M{fr) = lim ess sup f{u). 

iiS[r—<5,r-|-5] 


( 1 , 6 ) 
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Roughly speaking, F is obtained from / by “filling” the jumps of / in dis¬ 
continuity points. If / G where n > 1, the Filippov map 

F: R” —)■ 2®-" is defined as 

^(^) = n n COIN f{Bs{r)\N) (1.7) 

<5>0 m{N)=0 

where m is the Lebesgue measure and Bs{r) is the ball of centre r and 
radius 6. Of course F(ro) = f{ro) in all continuity points tq of /. Then 
to get existence in (1.3) one should replace / by F given by (1.7). If / is 
monotone and measurable then F is maximal monotone in R" x R” and 
locally bounded in R"', (see [1, Proposition 25]), and so, as shown in [2, 
Theorem 2.2], equation (1.3) has a unique strong solution (see also [3]). In 
the general case we consider here, the best that we can however expect is 
only a martingale solution for (1.3) (see Theorem 2.1, in which in general we 
do not have the uniqueness of the solution). 

The main existence result for equation (1.3) is established in Section 2, 
where it’s also given a “sliding mode” type result for this equation. 

In Sections 3, 4 and 5 it is studied a similar problem for equation (1.5) 
and also for a stochastic parabolic system. 

Notation We use the standard notation for the Sobolev spaces 
k = 1,2, Hq{0) and the Lebesgue integrable function spaces on O C R'^. 
The norm of Hq{0) is denoted by || • ||i and the norm of L^iO) by | • |p 
(1 < p < oo). The scalar product of L‘^{0) and the duality pairing between 
Hl{0) and the dual space H~^{0) is denoted by the same symbol (•, ■) 2 - We 
denote by F([0, T]; H) the space of all continuous 77-valued functions on [0, T] 
and we also refer to [4] for basic results pertaining stochastic processes with 
values in Hilbert spaces. Finally, we denote by F^(R), 7 = 0,1, the space 
of functions of class on R, with continuous and bounded derivatives up 
to order k. The norm in R or R"" is denoted by the same symbol | ■ |, the 
difference being clear from the context. 


2 Weak solution and “sliding” mode for the 
system ( 1 . 3 ) 

We shall study here system (1.3) where hF is a n-dimensional Wiener process, 
in a probability space (r2,F, P) and / G ^ Lip(R"'; £(R"; R"-)). 

We consider the Filippov map F: R"^ —)■ R” associated with / which was in¬ 
troduced in (1.7). 
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Definition 2.1. The system (12, P, (J^t)t>o, X)) is said to be a mar¬ 
tingale solution to (1.3) if (12, (J^t)t>o,P) is a filtered probability space on 

which it is defined an {iFt)t>o-Wiener process W and X is an {Xt)t>o-adapted, 
W^-valued, continuous process that satisfies P-a.s. the equation 

^it)+ j 'q{s)ds = x+ f cr(X(s)) (iiy(s), Vt > 0. (2.1) 

Jo Jo 

Here p G L°°((0,T) x 12), for each T > 0, is an {Xt)t>o-adapted process such 
that 

p G F{X), a.e. in (0, cxo) x 12. 

This definition extends verbatim to infinite dimensional equation (1.1). 
In literature such a solution is also called weak solution. A martingale 
solution which is Xfi' -adapted, where Xfi" is the completed natural filtration 
of W is called strong solution, see [5]. 

We have 


Theorem 2.1. Assume that f: R” —?■ R” is measurable and that 

|/(r)| < oi |r| -I- 02 , Vr G R"' 


( 2 . 2 ) 


where Oi > 0, 02 G R. Then for each x G R"" there is at least one martingale 
solution ipL, JF ,W, X) to (1.3) which satisfies the estimate 


E\X{t)\^ <Ct{\x\^ + 1), xeR^,te [0,T]. 

Proof. Consider the approximating equation 

dX, + ffiX,) dt = a{X,) dW, t G [0, T] 
A:(0) =x 

where /g is a smooth approximation of / given by 


(2,3) 


(2.4) 


(2,5) 


fe{'f')= / f ^ p{(^) dO , V£>0,rGR"'. 

Here p G (^“(R”) is any mollifier such that 

/ OO 

p{r) dr = 1. (2.6) 

•CX) 

Let Xg G L^(12; C'([0, T]; R"')) be the strong solution to (2.4). By (2.2) and 
Ito’s formula it follows that 


d \Xfit)\l < C \Xfit)\l dt + Xfit) ■ a(Xg(t)) dWt 


4 


and so by the Burkholder-Davis-Gundy theorem (see e.g., [4]) we have 


E 


sup \X^{t)\^ 
te[o,r] 


< G (1 + \x\‘^), 


> 0 


(2.7) 


(Here and everywhere in the following we shall denote by C several positive 
constants independent of e.) 

We set = (X^, W) and we consider = CiY^) (the law of Y^) that is 
VeiX) = G r] for each Borelian set F C G([0, T]; R") x G([0, T]; R”). Let 
us show that {ve} is tight in (G([0, T]; R”))^ = G([0, T]; R"') x G([0, T]; R"'). 
This means that for each <5 > 0 there is a compact subset F of (G([0, T]; R"'))^ 
such that z^£(F‘^) < <5 for all e > 0. We take for r > 0, 7 > 0, 

F = Br,, = {ye (G([0,T];R"))2 : \y{t)\ < r,Wt G [0,T], 

\yit)-y{s)\ <-f\t-s\2, Vf,sG[0,T]} 

Clearly, by the Ascoli-Arzela theorem, Rj ,,7 is compact in (G([ 0 , T]; R”))^. 
On the other hand, by (3.4) we have via Ito’s formula applied to the process 
t ^ \X,{t) - X,{s)\l 

I E|X,(t) - X,(s)|2 + E [\f,{X,{e)),XM - X,{s)) d9 

J s 

<CEj^\x,{e)\‘^de o<s<t<T. 

Taking into account estimate (2.7), we obtain via GronwalFs lemma that 

E\x,{t) - x,{s)\^ < c \x,{e)\^de < c \t - s\. ( 2 . 8 ) 

By estimates (2.7), (2.8) and by 


pF[\Y\> p] <E\Y\, Vp>0, 


we see that there are 7 , r independent of e such that as desired. 

Then by the Skorohod’s representation theorem there exist a probability 
space (G, X, P) and random variables X, IW, W such that C{X^, W^) = 
C{X^, H4) and for P-almost every ca G G 

W, ^ W, A, ^ A P-a.s. in G([0,T];R-) 
a(A,)^a(A) P-a.s. in G([0, T]; R^) ^ ^ 
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as e —)■ 0. We have also C{fs{Xe)) = C{fe{Xe)) and so by (2.2), (2.5) it 
follows that on a subsequence, denoted {Sn}, 

feni^er^) ^ V weak-star in L°°((0,T) X fi). (2-10) 

Let us show that 


rj G F{X), a.e. in (0,T) x 12. (2-11) 


We have by (2.5), 







e6) p{e) d9 G conYf{B^{X^{t,uj))), 

V(t, cu) G [0,T] X 12, 


and this implies that 

= { hm /,„(X,(2,a;))} C F{X{t,u)), W{t,uj) G [0,T] x 12. 

£n^0 

By (2.10), which implies of course also the weak convergence in L^((0, T) x 12), 
it follows by Mazur’s theorem (see e.g., [14, pag. 120]) that there is a convex 
combination of /e„, that is 


kfi 

i=l 

= 1, 0 < a'f'^ < 1, which is strongly convergent in L^((0,T) x 12) 
to fj and so on a subsequence again denoted {n}, 

lim ipn{t,u)) = fi{t,u)), a.e. it,u)) G (0,T) x 12. 

n—>-cx) 

Since hm„_j.oo <yi^n(t, cu) G F(Xe(2, cu)) we obtain (2.11) as claimed. 

If we dehne 


Ft = a{X,{s),W,{s)-,0 < s <t), 2 > 0, 

= (t(X(s), W(s); 0 < s < t), 2 > 0, 

then it follows that {We^Ft) and {W^F) are Wiener processes and that 
P-a.s., 

Xe{t)+ f UX,{s))ds = x+ f a{X,{s)dW,{s), V2G [0,T]. 

Jo Jo 
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Taking into account (2.10) and that P-a.s (see Lemma 3.1 in [8]) 


lim [ cr(X£(s)) (iiT£(s) 



a(X(s))dIT(s), 


Vte [0,T], 


we obtain that P-a.s. 


X(t) + 



X + 



a{X{s))dW{s), 


Wt e [0,T]. 


This means that the system ((1, P, hL(t),X(t)) is a martingale 

solution to (1.3). The estimate (2.3) follows by (2.7) which in turn implies 
that 

E|X,(t)|2 < ^(1 + l^n, Ve > 0, t e [0,T]. 

Such a process X can be extended to all of (0, cxo). □ 

Remark 2.1. If / and a are in L°°(R”) and 

n n 

(ffV)ij(x) >aJ2 il Vf = Ki) 6 E”, (2.12) 

ij=l 2=1 

then, as shown by A. Yu. Veretennikov [13], equation (1.3) has a unique 
strong solution X. (On these lines see also [8].) It should be said however 
that for the applications we have in mind, the nondegeneracy condition (2.12) 
is too restrictive. 


The sliding mode dynamics arises in differential systems with variable 
structure of the form (1.3) and a typical case is that when / has the form 


/(^) 


/i(r) if^(r)>0 

/2(r) if ^(r) < 0 


Vr G R 


(2.13) 


where g G 0^(R”), /i,/2 G 0^(R"',R"') and a G Lip(R"', £(R", R"')) satisfy 
the following conditions 


\fi(,r)\<aii\r\+a 2 i, Vr G R”, z = l,2 (2.14) 

sup {| V 5 '(r)| -I- |Z1^5((r)|} < cx) (2-15) 

reR" 

Xg{r) ■ /i(r) > a in {r G R"" : g{r) > 0} (2-16) 

V 5 '(r) • / 2 (r) < —a in {r G R"' : g{r) < 0} (2-17) 

\D^g{r)\\a{r)\‘^ <C*\g{r)\, Vr G R^ (2.18) 


where a > 0. We have 
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Theorem 2.2. Under assumptions (2.14)-(2.18) for each x G there is a 
martingale solution (11, P, IT, X) to (1.3) with the following properties: 
(i) if g{x) = 0 then P-a.s. g{X{t)) = 0, Wt > 0; 

(a) if g{x) 7 ^ 0 and r = infjt > 0 : g{X(t)) = 0} then 

P{t > t) <—{1 - e~^*)~^\g{x)\, Vt> 0 , (2.19) 

a 

where C = CiC*, Ci a positive constant independent of g and a. If C* = 0 
then 

P(r > t) < {at)~^\g{x)\, Vt > 0. 

Theorem 2.2 amounts to say that the manifold S = {x : g{x) = 0} is 
invariant for stochastic system (1.3) with / given by (2.13) and that for 
X ^ S the solution X have reached the manifold S by time t with a proba¬ 
bility greater or equal to 1 — {at)~^ \9{^)\- fhe classical automatic control 
terminology (see, e.g., [ 12 ]) this means that g{x) = 0 is a “sliding mode” 
equation for system (1.3) and S is a switching surface for this system. As 
a matter of fact this is typical “sliding” mode behaviour for the solution 
X = X{t) and its dynamics has two phases: the hrst phase is on time inter¬ 
val (0, r) until X reaches surface S and the second one for f > r in which 
X{t) evolves on sliding surface S. The reaching time r = t{u}) is a stopping 
time determined by (2.19). 

Proof of Theorem 2.2. We note hrst that the function / can be written as 

fir) = fiir)H{g{r)) + / 2 (r)hf(-^(r)), Vr G R” 

where H is the Heaviside function while the corresponding Filippov multi¬ 
valued function F (see (1.7)) is just 

F{y) = My)Higiy)) + My)Hi-giy)), Wy G R" 
and H is the multivalued Heaviside function 

{ 1 for r > 0 

[ 0 , 1 ] for r = 0 ( 2 . 20 ) 

0 for r < 0 . 

In the following we shall use the notations 


sgn(r) = H{r) — H{—r) 


1 

-1 


for r > 0 
for r < 0 


sgn(r) = H{r), Vr G R. 

Let X be the martingale solution to (1.3) given by (2.9) where / is as in 
(2.13). In order to prove the theorem we need a few apriori estimates on the 
solution to (2.4) which will be obtained by applying Ito’s formula to the 
function (j)\{u) = ipx{g{u)), where ipx G C^(R) is 


Va(0) = 0 

^xiy) = jy 
' ^xiy) = 1 + A 
^xiy) = -1 - A 


for I 2 /I < A, A > 0 
for y >2X 
for y < —2A 
for \y\ < 2A. 


( 2 . 21 ) 


On (0, cxo) such a function can be taken as 




r ^ 1/2 for 0 < // < A 
< (1 + A) 1 / for 2A < 1 / < cxo 
[pa( 2/) for a < // < 2A 


where px is a fourth order polynomial conveniently chosen and extend it by 
simmetry on (—cxo, 0). As a matter of fact, px is a smooth approximation of 
function y ^ \y\ and as easily seen 


I</^a(|/) - (sgn)A(l/)| < OA, V|/GR,|/^0, 


where (sgn)A is the Yosida approximation of sgn, that is 

A |r| for |r| < A 
(s'^)A(r’) = ■( 1 for r > A 

—1 for r < —A. 


We have therefore for all r G R \ {0} 

A^O 

Taking into account that, Vm,u G R"', one has 


lini (^(^(r) = sgn(r) 


Xpx{u) = p'x{g{u))X g{u), 

L>2(Pa(m)(u) = pl{g{u)){Xg{u) ■ v)Xg{u) + p'x{g{u))D‘^g{u){v) 


we obtain that 

dpx{g{X,m + p',{g{Xemfe{Xe{t)) ' Xg{X,{t)) dt 

=lTT[a%X,{t))a{X,{t))D^<j)x{Xemdt ( 2 . 22 ) 

+ aiX,it))dWit)-XMXs{t)), t G [0,T]. 
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Now, taking into account that on {g{Xi;{t)) ^ 0} 

= sgn{g{X,{t))) 

A^O 

and that in virtue of (2.18) and (2.21), 

ip'i{g{X,)) = 0 on{|(7 (X,)| >2A} 

\^'M^e))\<2 on{|^?(X,)|>2A} 

\D^g{X,)\\a{X,)\^ < C* |( 7 (X,)| in (0,T) x O 
\^'i{g{X,))\<C\-\ VA>0, 

letting A —)■ 0 in (2.22), we obtain that 

d\9{^e{t)) \ + ' V^(X£(t)) Sgn(^(X^(t))l[|g(x,(t))|>0] dt 

=i Tr[a*(X,(t))a(X,(t))D2^;(X,(t))]sgn(^7(X,(t)) dt (2.23) 

+ a{X,{t))dW{t)-Xg{X,{t)) 

By (2.5) and (2.13) we have 

UX,{t)) ■XgiX,it))sgnigiX,m 

j - 89) ■ vg{x,{t))p{e) de 

[g{Xe{t)-ee)>0] 

+ j f 2 {X,{t) - 89) ■ Vg{X,{t))p{9) d9 

[g{Xe{t)-ee)<0] 

and so taking into account (2.14)-(2.16) we get 

fe{X,{t)) ■ Vg{X,{t)) sgn{g{X,{t))) >a- (5(£)(1 + |X,(t)|) Vt > 0 
where (5(e) —)■ 0 as e —)■ 0. Taking into account (2.23) this yields 

d\g{X^{t))\ + a l[| 3 (x,(t))|>o](l - S{8)\X^{t)\) dt 

<i TT[a*{X^{t))a{X^{t))D‘^g{X^{t))]sgn{g{X^{t))) dt 
+ lMxpm>o]\9{Xem-^Xg{X,{t)) ■ a{X,{t)) dW{t) 
and therefore, for 0 < s < f < T, we have P-a.s. 

|9(A'.(«))I +a ^‘ll[wx.«,))|>0l(l - ^(<r)|A'A»)l)<i» 

<l9(AAi>))|+ f\MxAm>»MMO))\-Xa(XA0)) ■ a(XM)dW(e) 

J s 

+ I j\T[a*{XM)^{Xe{0))D^g{X,{9))HnmXe{0)))d9. 
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The same inequality remains of course true for [X^, W) and so letting e —)■ 0 
we get that for X given by (2.9), we have 


[ ^[\g{x(e))\>o]^^ 

J S 

<|9(V(i,))| + l^‘Tl'[(r-(X(«))(r(X(«))C 29 (X(»))]sgn(( 9 (X(»)))<i« 

+ / l[| 5 (x(e))|>o]V^(^W) • ^(^(^)) dW{e), 0 < s < f < cx>, P-a.s. 

J S 

Taking into account (2.18) we get 


| 5 '(^(^))I+® f '^[\g{x{e))\>o]dd 

J S 

+ [‘ 1 -|I»(X(»))|> 0 | Vs(V(«)). a{xm dW{f)), 

J s 

and so by the Gronwall lemma 


<C |^/(X( 0 ))| 


e + O' j e ^^l[\g(^x{d))\>o]dd < e ^^| 5 f(X(s))| + 

(2.24) 

^~^%\,ixm\>o]^9ix{e)) ■ aix{e)) dw{e), o < s < t < oo. 

In particular, it follows by (2.24) that if g{x) = 0 then g{X{t)) = 0 P-a.s. 
for all f > 0. Moreover, by (2.24) it follows that Z{t) = \g{X{t))\e~^^ is a 
nonnegative super-martingale and therefore for any couple of stopping times 
Ti < T 2 we have Ziji) > Z{t 2 ). This implies that if r = inf{t > 0 : \Z{t)\ = 
0} we have that Z{t) = Z{t), P-a.s. for t > r. On the other hand, by (2.24) 
and (2.18) it follows that 



EZ(t) + a e-^^P{T>s)ds<\g{x)\+C EZ{s)ds, Vt > 0 

Jo Jo 

and therefore 


C 


P(r > f) < — (1 - e ^ ) |^(a:)|, Vf > 0, 

a 


which is just (2.19). This shows that X{t) reaches the manifold S in stopping 
time r and remains there for t > r with a probability P greater or equal 
^ (1 — |5'(^)I- The proof is complete. □ 
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Remark 2.2. If conditions (2.16), (2.17) are satisfied with a = 0 in Theorem 
2.2 then only part (i) follows. 

Theorem 2.2 can be used to design feedback controllers for stochastic 
differential systems with a sliding mode dynamics on a given surface S = 
{x : 5 f(a:) = 0}. Such an example is presented below. 

Example 2.3. Consider the controlled stochastic second order system 

X + Oi X = (To(X, X)/3 + M in(0,oo). (2.25) 

We assume that ao € Lip(R^). 

Our aim is to hnd a feedback controller u = —fo{X,X) such that the 
corresponding closed loop system 

X + ai X + /o(X, X) = ao(X, X)$ 

X(0) = xo, X(0) = xi ^ ^ 

has the sliding mode equation 

a2X + X = 0, (2.27) 

for some 02 G R. Here /3 is a Brownian motion in a probability space (O, X, P) 
and $ is the associated white noise. 

We choose 


/o(nX2) = asgn(a2ri Trs), V(ri,r2)GR^ (2.28) 

where a > 0 and rewrite equation (2.26) as 

dXi — X 2 df = 0 , . 

<7X2 + di X2 dt a sgn(a 2 Xi + X2) dt = (To(Xi, X2) dj 3 

for t > 0 , where as us usually sgn(M) = for u 7 ^ 0 . 

Equation (2.29) is a “jump” system of the form (1.3) where 

= ( a, r. + «sgu(a.r. + r,) ) ’ ^ 

and so / is of the form (2.13) where 

/>M=(a,7 + „)' = ( a. rT-a ) ■ = (’’i.^ 

g{r) = a2ri+r2, r = {ri,r2). 

It is easily seen that conditions (2.14)~(2.18) hold and so Theorem 2.2 is 
applicable to the present case. We get 
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Corollary 2.3. The stochastic closed loop system (2.29), equivalently (2.26), 
(2.28), has the “sliding mode” (2.27). More precisely, for every (a;o,Xi) G 
there is a martingale solution {Xi(t), X 2 {t)) which reaches the surface S = 
{(xi, X 2 ) : a 2 Xi + X 2 = 0 } in time t with a probability > 1 — (a t)“^|a 2 Xo + Xi|, 
and remains P-a.s. on this surface after that time. 

This describes a typical “sliding-mode” behaviour for solutions X to 
(2.26), namely 

aiX{t) + X{t) = 0 

on (to, cxo) X Qq where P{Tlo) > 1 — {atQ)~^\a 2 XQ + xi\. (We refer to [9], [10], 
[ 11 ], for references and other significant results on “sliding-mode” behaviour 
of stochastic differential systems). 


3 Existence of a weak solution to heat equa¬ 
tion ( 1 . 4 ) 

The following hypotheses will be assumed throughout in the sequel. 

i) / G L[^JR) and |/(r)| < ai|r| -|- hi, Vr G R 

ii) W is the cylindrical Wiener process (1.2) where is an orthonormal 

basis in L‘^{VL) given by — Ae^ = XjCj in O] Cj = 0 on dO and 

CX) 

(3.1) 

i=i 


iii) b G C^(R)nLip(R). 

Definition 3.1. Let x G L‘^{0). We call weak (martingale) solution to 
( 1 . 1 ) a tuple {VL,X,{Xt)t>o,^,W,X), where ( 12 , (J^i)t>o, P) is a filtered 

probability space where there are defined a {Xt)t>o-Wiener process W and a 
continuous {Xt)t>o-adapted L‘^{0)-valued process X = {X{t))t>o such that, 
P-a.s., 

X{t) = e-^^x+ [ e-(*-*)^r/(s) ds + [ 6(X(s)) diy(s), (3.2) 

Jo Jo 

where rj G L°°((0,T) x O x Q) is a {Xt)t>o-adapted process such that 

f] G F{X), a.e in (0,T) x O x kl. (3.3) 
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Here A = —A with V(A) = 11^(0)1111^(0) and is the Co-semigroup 
on L‘^(0) generated by —A. 

We note that the linear operator h(X) arising in (3.2) is dehned by 

OO 

h(X)h = J2f^jKX){h,ej)2ej, Mh e L‘^(0). 

i=i 

The construction of a weak (martingale) solution. We consider the 
approximating equation 

dX, - XX, dt + f,(X,) dt = b(X,) dW, in (0, T) x O 
X, = 0, on(0,T)xdO 

x,(0,O = x(O, eec 

where £ > 0 and, as in the hnite dimensional case (see (2.5)), 


(3.4) 


fe(r) = - f{s)p{!-^)ds = f(r-Eff)p(e)de, Vr 6 E, (3,5) 


Here p E C“(R) is such that 

p(o) > 0 p(e) = p(-9), p(e) = 0 for \e\ > i. 


p(e)d9 = l. (3.6) 


Clearly by (i) we have 

/e e C^(R), \fe(r)\ <ai\r\+bi + aie, VreR,e>0. (3.7) 

By standard existence theory for inhnite dimensional stochastic equations 
with Lipschitz nonlinearity it follows that (3.4) has a unique strong solution 

A, e L2(H; C([0, T]; L\0))) 1 ^^(0, T; HI(0))), (3.8) 

see [4, pag.45]. By Ito’s formula we get P-a.s. 


\\X,(t)\l+ f\\X,(s)\\lds+ [\f,(X,(s)),X,(s)) 2 ds 

'o Jo 

ft °° rt 

^ I . / ■ I _ — t— — I ^ I/. I/ll .A 

'0 


Pt 

+ l P^j\K^e(s)) 

Jo j=l 


ej\lds+ I {b(X,(s))dW(s),X,(s)) 2 ds, 

vte [0,T], 


and so by the Burkholder-Davis-Gundy formula we obtain by some calcula¬ 
tion involving (i)~(iii) 

E sup \X,(t)\l + 'E f ||Ae(s)||^ ds < C (|a:| 2 -I-1), Ve > 0, (3.9) 

te[o,T] Jo 
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where C is independent of e. By (3.7) we also have 


E snp \f^{X^(t))\l<C {\x\l + l). 

iG[0,T] 


Then on a snbsequence, again denoted in the same way, we have for £ —)■ 0 

(3.10) 


^ X weak-star in L'^(0, T; L‘^{0)) 
weakly in L^{Q,T-Hl{0)) 


and 


fe^rj weak-star in L°°((0, T); (hi; L^(C>))) (3-11) 

X,{t) X{t) weakly in ^^((0, T) x O)) (3.12) 

h{X,) ^ b* weakly in L‘^{n; L‘^{{0, T) x O)) 

dX — AX dt + rjdt = b* dW in (0, T) x O, 

X(0)=a; in C>, (3.13) 

X = 0 on (0,T) X (90, 


that is 


X{t) — / AX{s)ds+ / ri{s) ds = X + / b*{s)dW{s), VfG[0,T], P-a.s. 

Jo Jo Jo 

Assume now that x G Hq{0). Then by an application of Ito’s formula in 
(3.4) to the function a; —| ||x||f we get for some Ci, C 2 > 0, 

l\\XM\\l+ f\AX,{s)\lds<l\\x\\l + C, [\\AXM\l + \X,{s)\l)ds 
Jo Jo 

pt ^ pt 

+i / Tf^^jMX,)ej\\lds+ / {Ab{X,{s)),dW{s)) + C 2 
Jo Jo 

and in virtue of (ii), (iii) this yields via Burkholder-Davis-Gundy formula 

E sup ||A:,(t)||2+ E / |AX,(s)|2ds < 0(11x11? + 1), Ve > 0. (3.14) 

te[o,r] Jo 


(Everywhere in the sequel we shall denote by C several constants independent 
of e.) 

Then, in this case besides (3.10), (3.12), we also have 

X G C{[0,T]] L\0)) n L°°{0,T] H^{0))) n L\0,T] H\0))). 

(3.15) 
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Since the weak convergences (3. 10 )-(3. 12 ) are not sufficient to conclude that 
(3.3) holds, then proceeding as in the proof of Theorem 2.1 we shall re¬ 
place {Xg} by a sequence of processes dehned in a probability space 

{fl, P, IT} such that C,{X^) = C{X^) where C is the law of the process. 

To this end, consider the sequence {i^e}e>o of probability measures, = 
C{Xe), that is = P(X£ G B) for any Borelian set B C C([0, T]; L‘^(0)). 

We have 

Lemma 3.1. Let x G HliO). Then the sequence {i'e}e>ei is tight in the space 

C(|0,T];n(O)), 

Proof. This means that for each 5 > 0 there is a compact subset B of 
C{[0,T]; L'^{0)) such that UeiB^) < h for all £ > 0. We take for r > 0, 
7 > 0, 

B = Br,, = {ye C{[0,T]-L\O)) : \y{t )\2 < r,Wt G [0,T], 

\\y\\L^(o,T-,Hpo)) < r, \y{t)-y{s )\2 <-f\t-s\2, yt,s G [0,T]} 

Clearly, by Ascoli-Arzela theorem, Br^^ is compact in C([ 0 , T]; L^(0)). On 
the other hand, by (3.4) we have via Ito’s formula applied to the process 
t ^ \X,{t) - X,{s)\l 

\ E|x,(f) - x,(s)|2 + E J\vx,{e), V(X,( 0 ) - x,(s )))2 de 

+Ej^{f,{x,{e)),x,{e)-x,{s)) 2 de<CE\x,{e)\lde o<s<t<T. 

Taking into account estimates (3.7), (3.9) we obtain via Gronwall’s lemma 
that 

E|A,(t) - x,(s)|2 < c j\\x,{e)\l + \vx,{e)\l) de<c\t-s\. ( 3 . 16 ) 

By estimates (3.9), (3.14), (3.16) and taking into account that 

pE[\Y\> p] <E\Y\, Vp>0, 

we infer that there are 7 , r independent of e such that < 6, as 

desired. □ 

Then by the Skorohod theorem (see, e.g.. Theorem 2.4 in [4]) there are 
a probability space (0,J^, P) and the stochastic processes X, {Xs}e>o on 
(O, P) such that the law C{X^) of W coincides with C{Xs) and P-a.s. 

mC{[0,T]-L\O)) (3.17) 
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as e ^ 0. We have also C{X) = C{X). Since C{fe{Xi;)) = ^(/(.(Xg)), 
C{a{X^)) = C{a{Xs)) by (3.17) and (3.5) we see that 


b{X,)^b{X), 


a.e. in (0,T) x O x Cl, 


where £(?)) = C{r]) and 

fj G F{X), a.e. in (0,T) x (9 x 11. 


(3.18) 


(3.19) 


The latter follows as in the proof of Theorem 2.1 taking into account that in 
this case F is given by (1.5), but we omit the details. 

We set 

M,{t) = X,{t)-x - [ AX,{s)ds + 

Jo 

and 

M,{t)=X,{t)-x- [ AX,{s)ds+ [ fs{Xe{s))ds, t e [0,T]. (3.21) 

Jo Jo 

It turns out that is a square integrable martingale on (f2,X, P) with 
respect to the hltration Ft = clXs; s < t} because since C{M^) = C{M^) 
and Mg is a square integrable martingale on (fl, F, P) we have 


MX,{s))ds, te[0,T] (3.20) 


E[{X,{t) - X,{s)- 


Ax,{e) d 9 +MMO)) de)x{x,)] = o (3.22) 


for any bounded continuous function x and all 0 < s < f < T. 

Passing to the limit in (3.20) and taking into account (3.17)-(3.19) one ob¬ 
tains that the process 


M{t) =X{t)-x- 


AX{s)ds+ / fi{s)ds, 


t > 0 


is an L^(C>)-valued martingale with respect to hltration Ft = cr{X(s), s <t}, 
t G [0,T], with hnite quadratic variation, see [4, pag.234]. Then by the 

representation theorem 8.2 in [4] there is a larger probability space (12, F, P), 
a hltration {Ft\t>o and an L^(C>)-cylindrical Wiener process W{t) on it such 
that, P-a.s., 

M{t)= f b{X{s))dWit), fG[0,T]. 

Jo 

This means that the system (12, X, {Xt}t>o, P, lP(t), X(f)) is a martingale 
solution to (1.1). We have proved therefore 
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Theorem 3.2. Under Hypotheses (i), (ii), for each x G Hl{0), there is at 

least one martingale solution (fi, P, X) to equation (1.1) and X 

is given by (3.17). Moreover, we have 

X G L°°(0, T; H^{0))) n L\0, T; H\0))). (3.23) 

We note that (3.23) follow by (3.14) and (3.17). 

Remark 3.1. Under additional assumptions on b (for instance if it is in¬ 
dependent of X) it turns out that the martingale solution X is the unique 
strong solution, see [5]. (See also Remark 2.1) 


4 Sliding mode control of the stochastic heat 
equation 


For parabolic stochastic equations of the form (1.1) a “sliding” mode dynamic 
arises for discontinuous (“jump”) functions / : R —?■ R of the form (2.13), that 

is 


f{r) 


/i(r) for 5 ((r )>0 

/ 2 (r)) for g{r) < 0 ’ 


(4.1) 


where g, /i, /2 are given continuous functions. 

As in the previous hnite dimensional case, the objective of the “sliding- 
mode” control is to design for the linear time invariant system 


dX — AX dt = du in (0, T) x O 
X = 0 on{0,T)xdO (4.2) 

X{0) = X in O 


a stochastic feedback controller of the form 


du = -f{X)dt + b{X)dW (4.3) 

such that the “sliding” motion occurs on the manifold S = {X : g{X) = 0} 
which is also referred as “sliding” or “switching” surface. Roughly speaking, 
this means that any trajectory of the closed loop system (4.2)-(4.3), which 
starts from initial state x, reaches the sliding surface S at a certain time 
to and remains there for t > to- As a matter of fact, this last phase of the 
dynamics is called “sliding mode”. 
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Of course in virtue of Theorem 3.2 a weak solution X to (4.2) in the sense 
of Dehnition 3.1 exists for the extended multivalued closed loop system 


dX -AX dt + fi{X)H{g{X)) dt 

+f 2 {X)H{-g{X)) dt = a{X) dW 

X = 0 

X(0) = X 


in (0,T) X O 
on (0, T) X dO 
in O. 


(4.4) 


Here H is the multivalued Heaviside function (2.20) on R. To begin with we 
shall prove hrst an invariance result for the manifold S = {X : g{X) = 0}. 

Theorem 4.1. Let g G 0^(R), fi, /2 be continuous functions which satisfy 


assumption (i) and let b, satisfying (hi), be such that 

h\r){gg’' + {g'f){T) < Cg^r), Vr G R (4.5) 

g{r) g'\r) + {g'{r)Y >0 Vr G R (4.6) 

/i(r) g\r) > 0 for g{r) > 0 (4.7) 

/ 2 (r) g'{r) < 0 for g{r) < 0. (4.8) 

for some C > 0. Then, for all x G Hl{0) such that g{x) = 0 on O, there is 

a martingale solution (0,7^, P, X) to system (4.4) such that 

g{X{t)) = 0, Vf G [0,T], P-a.s.. (4.9) 

Proof. We start with the approximating equation (3.4). We apply the ltd 
formula to function x g‘^{x) and get 


d / g\X,it,0)df+ 2 / {gg'' + {gr)iX,it,0)\^X,it,Ordfdt 


lo 


'o 


+ 2 / MX,it,0)9iX,{t,0)9\X,it,0)dfdt = 
Jo 

OO « 


i=i 


'o 


+ I b{X,{t,0)9{X,{t,0)9'{X,{t,0)e,dfdfJ,{t) 

i=i 


'O 
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Taking into account (3.5),(3.7), we obtain that 

I f,{XMx,)a'(x,)d( = j p{0) 

( j h(X,-Ee){gg'){X,-Ee)df+ 

[g{XE-ee)>0] 

j UX,-e9){gg'){X,-e9)d^')d9 + Ut), Vt G [0,T], 

[g{Xe-ee)<0] 

where 

Cs{t)<Ce [ {\X,{t,0\ + l)d^ 

Jo 

with C = CiC^ Cl > 0; thus, by (4.5)-(4.7), this yields 

E [ g\X,{t,0)d^<CE [ [ g\X,{t,0)d^ds 
Jo Jo Jo 

+ 5{e)E [ [ {X,{t,0^ + l)d^ds, VtG[0,T]. 

Jo Jo 

where lim£^o<5(£) = 0 and the constant C is independent of £. 

This yields via Gronwall’s lemma, 

E [ g\X,{t,0)d^ < 6{e) exp(C't), Vt G [0,T]. (4.10) 

Jo 

If Xg is dehned as in the proof of Theorem 3.1, that is >C(Xe) = C{X^) and 
(3.17) holds, we get by (4.10) that 

E [ g‘^{X,{t,i))di < 5{e) exp(C'f), Vt G [0,T], Ve > 0 
Jo 

where E is the expectation in probability space (fi,X, P). Hence, letting 

e tend to zero we get g‘^{X) = 0, dt x d^ x P-a.e. in (0,T) x O x Q as 
claimed. □ 

Remark 4.1. In the particular case where the function 

F(r) = /i(r)H(^(r)) + f 2 {r)H{-g{r)), r G R 

is a maximal monotone graph in R x R with R{F) = D{F) = R, equation 
(4.4) has a unique strong solution X. This happens for instance if /,, i = 1, 2, 
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are monotonically nondecreasing continuous functions such that fi > f 2 on 
R. and g{x) = x (see [2]). Then the corresponding system (4.4) 

dX -AXdt+ {fi{X)H{X) + f2{X)H{-X)) dt = 

b{X)dW, in(0,T)xC> (4.11) 
X = 0, on(0,T)xaC>, 

has the invariant manifold X = 0 . 

Under stronger assumptions on g and b it turns out that the closed loop 
system (4.2)-(4.3) (equivalently (4.11)) has a “sliding” mode dynamics with 
the switching manifold S = {X : g{X) = 0}. Namely, we assume that 


Hypothesis 4.2. fi, i = 1,2 satisfy assumption (i) and g E U^(R), b G 
C'^(R) n LipfR) are such that 

\g'\r)\ |6(r)|2 < C* \g{r)\, Mr e R, (4.12) 

fi{r)g\r) > a if g{r) > 0; f 2 {r) gfr) < -a if g{r) < 0 (4.13) 

g',g"EL^{R), g"sgnig)>0 on R, (4.14) 

where a > 0 . 

We note that by Theorem 3.2, equation (4.11) has a martingale solution 
X given by (3.17). 

Theorem 4.3. Under Hypothesis 4.2, for each x G Hl{0) there is a mar¬ 
tingale solution (12, X, P, hU, X) to (4.2), (4.3) such that for r = inf{t : 
\g{X{t))\ = 0} we have 

P(r >t) < —(1 - e~^^)~^ \g{x)\^- (4.15) 

a 

where C = CiC*. Moreover, if g{x) = 0 a.e. in O, then g{X(t)) = 0 for all 
t > 0. 

Proof. The proof is very similar to that of Theorem 2.2, so it will be sketched 
only. If Xe is the solution to equation (3.4) and Xg such that C{X^) = C{X^), 
we get via Ito’s formula applied to function x —)■ ip\{g{x)) and after letting 
A ^ 0 

d |^(X^(t ))|2 + < {b{Xe{t))dW(t), g'{X,{t)) sgn{g{X,{t )))2 

OO „ 

fcl Jo 
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where ip\ is the function introduced in the proof of Theorem 2.2 (see (2.21)). 
We have used here condition (4.14) which in virtue of (2.21) yields 




A^O 


(4.16) 

= -lini / VX,(t,O)-V(^',(r?(X,(f,O))^ 7 '(X,(f,O))c^e< 0 , 


o 


and also assumption (4.13) which, as we have seen in the proof of Theorem 
2 . 2 , implies that 

UX,)g\X,)sgng{X,) >a- <5(£)(1 + |X,|). 

Now using (4.12) and letting e —)■ 0 we get 

| 5'(-^(^))|2 + a ^ l[|g(X(6»))|>0]^^ - Ifi '(^(^)) l 2 + | 5'(-^(^))|2 d 6 *+ 

[\biX{e))dW{e),g'{X{e))sgn{g{X{e )))2 for 0 < s < f 


for some constant C > 0 . 

This yields (see (2.24)) 

e-^^\g{X{t))\2 + a j < e-^^\g{X{s))\2+ 

J S 

J s 

0 < s < t < oo. 


Here Z = \g[X[t))\ 2 e is a nonnegative supermartingale and so Z{t) = 
Z(t) P-a.s. for f > r = inf{f > 0 : \Z{t)\ =0}. Taking expectation, we get 

[ FjZ{s)ds 
0 

which implies the desired estimate (4.15). □ 

By (4.14) we see that g is convex on [^f > 0], concave on [g < 0] and that 
[r : g{r) = 0] = [ai, a 2 ] is a closed interval. 

Hence the switching manifold S = {X : g{X) = 0} is of the form 

E = {X G L'^iO) : ai < X < as} 

and so under the assumptions (4.12)-(4.14) the closed loop system (4.2)- 
(4.3) has for each x G Hl{0) a martingale solution which reaches the set S 
(that is the interval [ 01 , 02 ]) by time t with probability estimated by (4.15) 
and remains in this interval after that time. 


EZ(f) + o / e ‘^®P[r > s]ds < \g{x )\2 + C 
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5 “Sliding” mode control of a stochastic pa¬ 
rabolic systems 


Consider here the parabolic system 

dX - AX dt + /i(X, Y) dt = bi{X, Y) dWi, 
dY - AY dt + f2{X, Y) dt = b2{X, Y) dW2, 

X(0) = a:(O, b^(0) = |/(O, 

X = Y = 0, 

where /j G C'(R^) satisfy assumption (i), bi G C'^(R) fl Lip(R), i = 1,2 and 
Wi, W 2 are Wiener processes of the form (1.2) in the space H = Li^iO) x 
L'^{0). Let R2 ^ R, ^ G ^^(R^) be given. 

Arguing as in the proof of Theorem 3.2 it follows that for each {x,y) G 
Hh{0) X HliO), system (5.1) has a martingale solution {X,Y) obtained as 
limit of solutions {X^, Y^) to corresponding approximating system 


in (0,T) X O 
in (0,T) X O 
e G C> 

on (0,T) X dO 


d 


X 

Y 


+ 


-A 

0 






dt 




fbi{X,Y)dWi\ 

\b2{X,Y)dW2) 


where F: R^ ^ R^ is given by 


(5.2) 


F 


fl in [g > 0], 
/2 in [g < 0]. 


and 

Fe{ri,r2) = 
Assume further that 



p(r — e9) F{6) dO, 


(ri,r2) G R^, e > 0. 


\D‘^g{r)\{\bi{r)\^+ \b 2 {r)\^) <C\g{r)\, Vr G R^ (5.3) 

/i(’^)5'ri(’^) > a in [r : 5 ((r) > 0] (5.4) 

f 2 ir) gr 2 ir) < -a in [r : g{r) < 0] (5.5) 

where a > 0 and ((?ri, 5 'r 2 ) = '^g, 

gr,r, Sgng > 0, - ^nn^rara) Sgn^ < 0. (5.6) 


We have 
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Theorem 5.1. Under assumptions (5.3)-(5.6) for each {x,y) G Hq{0) x 

Hq{ 0) there is a martingale solution {12, iF, P, W, {X, F)} to (5.1) such that 
if T is the stopping time r = infjt : g{X{t),Y{t)) = 0} then 

]P[r >t]< ^(1 - e“^*)"^|^(x,|/)|(L2(o))2 (5.7) 

for some constant C > 0. 

The proof is exactly the same as that of Theorem 2.2 where the approxi¬ 
mating equation (3.4) is replaced by (5.2). We note that in this case the 
corresponding inequality (4.16) is a consequence of hypothesis (5.6). The 
details are omitted. 

A particular example is 

giri,r2) = apri + a2r2, Vri,r2eR 

which, for fi and /2 satisfying condition 

aifi{r)>a in {oi ri 02 r 2 > 0} 
a 2 f 2 {r)<-a in {oi ri-h 02 r 2 < 0} 

where a > 0, imply that system (5.1) has a martingale solution {X,Y) that 
reaches the linear manifold 

S = {«! X + a2Y = 0} 

in a time t with probability P > 1 — Ct~^\aiX + 0^2 y\{L'2{o)P and remains 
on this manifold afterwards. 
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